Direct observation of charge order and orbital glass state in multiferroic LuFe 2 04. 



o 
o 



o3 



i 

c 

o 
o 



> 
O 

m 

oo 
o 



A.M. Mulders 1,2 , S.M. Lawrence 1 , U. Staub 3 , M. Garcia-Fernandez 3 , 

V. Scagnoli 4 , C. Mazzoli 4 , E. Pomjakushina 5 * 6 , K. Condor 5 , Y. Wang 7 . 
1 Department of Imaging and Applied Physics, Curtin University of Technology, Perth, WA 6845, Australia 
2 The Bragg Institute, Australian Nuclear Science and Technology Organization, Lucas Heights, NSW 2234, Australia 
3 Swiss Light Source, Paul Scherrer Institut, 5232 Villigen PSI, Switzerland 
4 European Synchrotron Radiation Facility, BP 220, 38043 Grenoble Cedex 9, France 
5 Laboratory for Developments and Methods, Paul Scherrer Institut, 5232 Villigen PSI, Switzerland 
6 Laboratory for Neutron Scattering, ETHZ & PSI, 5232 Villigen, Switzerland and 
7 Department of Applied Physics, The Hong Kong Polytechnic University, Hong Kong, China 

(Dated: July 17, 2009) 

Geometrical frustration of the Fe ions in LuFe204 leads to intricate charge and magnetic order and 
a strong magnetoelectric coupling. Using resonant x-ray diffraction at the Fe K edge, the anomalous 
scattering factors of both Fe sites are deduced from the (h/3 fe/3 1/2) reflections. The chemical shift 
between the two types of Fe ions equals 4.0(1) eV corresponding to full charge separation into Fe 2+ 
and Fe 3+ . Polarization and azimuthal angle dependence of the superlattice reflections demonstrates 
the absence of differences in anisotropic scattering revealing random orientations of the Fe 2+ orbitals 
characteristic of an orbital glass state. 

PACS numbers: 71.70.Ch; 75.40.Cx; 78.70.Ck 



New materials that exhibit strong magnetoelectric cou- 
pling are fascinating because a large coupling between 
ferroelectric and magnetic interactions is rare, and its ori- 
gin often unclear. Competing interactions lead to novel 
ground states that give rise to unusual material proper- 
ties, i.e coexistence of spontaneous magnetic and ferro- 
electric order [lj]. The ability to control electric polar- 
ization with a magnetic field or the magnetization with 
an electric field [2, H, H| makes these multiferroic ma- 
terials promising candidates for novel applications such 
as 4-state memory and switchable magneto-optical de- 
vices. In most ferroelectric materials, electric polariza- 
tion arises from covalent bonding between anions and 
cations or the orbital hybridization of electrons. Alterna- 
tively ferroelectric polarization may arise from frustrated 
charge order as reported for LuFe204 This compound 
is of particular interest because, in addition to ferroelec- 
tricity, magnetism originates from the same Fe ions and 
this holds the promise of strong magnetoelectric coupling. 
The ferroelectric and magnetic order take place at and 
near ambient temperature which provides the potential 
for room temperature multiferroics. 

The crystal structure of LuFe2C>4 consists of a triangu- 
lar double layer of iron ions, forming trigonal bipyramids 
with five oxygen nearest neighbors, in which an equal 
amount of Fe 2+ and Fe 3+ arc believed to coexist at the 
same site Q. The occurrence of different charge order 
schemes has already been studied in more detail theoret- 
ically and experimentally 0,0]- LuFe204 adopts a ferro- 
electric ground state below ~350 K, while below ~250 K 
two-dimensional magnetic order is established in the tri- 
angular planes which enhances the ferroelectric polariza- 
tion by 20%, illustrating coupling between the magnetic 
and ferroelectric order The observation of the (1/3 

1/3 13/2) reflection using resonant x-ray Bragg diffrac- 



tion (RXD) further supports the existence of charge or- 
dering [Icj ]. 

RXD has become a powerful technique to study charge, 
orbital and magnetic arrangements. Tuning the energy 
of the incoming radiation to an absorption edge permits 
recording Bragg reflections with enhanced sensitivity to 
the specific ion and its electronic configuration. In case 
of the Fe K edge, the incident x-rays virtually excite an 
electron from the Is core level to the empty 4p states, 
followed by a decay of the electron back to the core hole. 
This effect results in a significant variation in the atomic 
scattering factors of the Fe ions for x-ray energies close 
to the Fe K edge. The atomic scattering factors are also 
affected by variations in charge state, ordered aspherical 
electron densities or ordered magnetic moments. In par- 
ticular asphcricity of the atomic electron density results 
in anisotropy of the tensor of x-ray susceptibility (ATS). 
Each of these phenomena has a specific dependence on 
the polarization of the incoming radiation and the orien- 
tation of the sample with respect to the scattering geom- 
etry The significance of RXD has been demonstrated, 
among others, in the manganites [HI], nickelates 12 1 and 
magnetite [13| . 

In a perfectly charge ordered state, each site may be 
considered as having an excess and a deficiency of half an 
electron respectively, compared to the average ion valence 
of Fe 2 5+ . Alternatively charge disproportionation with 
fractional charges may exist as exemplified in nickelates 
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151 ]. Frustration arises because every excess charge 
prefers a deficiency charge as a neighbor which is not 
possible on a triangular lattice. However, in the presence 
of a second triangular layer a net transfer of charge from 
the first layer to second occurs because then it is possible 
to have charge order on each layer in a honeycomb lattice 
arrangement HI, 17j ]. 
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The crystal field of the trigonal bipyramids splits the 
3d states of LuFe204 into two doublets (d xy /d x 2_ y 2 and 
dxz/dy Z ) and a singlet (d z 2) [jjfr . Fe 3+ with five 3d elec- 
trons is spherical while Fe 2+ with six 3d electrons ex- 
hibits doubly degenerate orbital degree of freedom in the 
d xy /d x 2_ y 2 ground state. 

In this paper we present RXD data with azimuthal an- 
gle and polarization analysis and our results clarify that 
the charge order is close to electronic states of Fc 2+ and 
Fe 3+ , in contrast to small values of disproportionation 
observed in nickelates [lj| and manganites 11 1. More- 
over the absence of scattering due to ATS demonstrates 
a glass state of the Fe 2+ orbitals in agreement with cal- 
culations (ljj ]. 

Polycrystalline LuFe2C>4 was prepared by a solid state 
reaction as reported in ref. [1 91 ] . Starting materials 
of L112O3 and Fe2C>3 with 99.99% purity were mixed, 
pressed into pellets and sintered at 1200°C during 6 
h in H2/He/C02 atmosphere (H2/CO2 ratio 1/3) and 
quenched into ice water. After grinding, the obtained 
powder was hydrostatically pressed and sintered at the 
same conditions during 3h. The crystal growth was car- 
ried out using Optical Floating Zone Furnace with four 
1000W halogen lamps as a heat source, growth rate 1 
mm/h, 2 bar pressure of CO2/CO mixture (5/2 ratio). 
The single crystal was cut and samples have been pol- 
ished perpendicular to the [001] and [110] directions. The 
magnetic ordering temperature was determined with a 
SQUID magnetometer at 240 K and pyroelectric current 
measurement confirmed ferroelectric order below ^330 K 
and enhanced ferroelectric order below ~220 K. 

Various (h/3 k/3 1/2) superlattice reflections were 
recorded at the Fe K edge at beamline ID20 of the ESRF 
[20I ]. Polarization analysis was performed using a MgO (2 
2 2) analyser crystal of which the polarization efficiency 
was determined at 0.98 for the energy of the Fe K edge. 
In addition, RXD was recorded at the MS beamline at the 
SLS (see also supplement 2l|) using the Pilatus 2D de- 
tector [22| . The background, mainly originating from the 
fluorescence of the sample, was determined from selected 
border regions of the area detector and subtracted. The 
integrated diffracted intensity was corrected for polariza- 
tion, absorption and sample geometry. The absorption, 
n(E), was obtained from the (006) reflection by itera- 
tion of the calculated anomalous intensity according to 
spacegroup R3m and fx(E) deduced from the ratio be- 
tween calculated and integrated intensity recorded with 
the 2D detector. Besides n(E) was recorded at the X- 
ray Absorption Spectroscopy beamline at the Australian 
Synchrotron using powdered LuFe2C>4 pressed with cel- 
lulose and Fe foil as energy calibration. It confirmed the 
validity of the method to obtain n(E) from the (006). 
RXD recorded at the SLS was used to analyse the energy 
dependence while RXD recorded at ID20 was used to in- 
vestigate polarization dependence and rotation about the 
scattering vector q (azimuthal angle iff). "J is defined zero 
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FIG. 1: Integrated resonant diffraction intensity of the charge 
order (2/3 -1/3 1/2) reflections recorded with 2D detector at 
T=15 K (see inset top left for I — 11), corrected for back- 
ground, polarization, absorption and sample geometry, nor- 
malized (shifted for clarity) and compared to charge order 
model described in the text. The inset (top right) shows that 
the superlattice intensity gradually disappears at the ferro- 
electric ordering temperature. For A, B, C and D see text. 



when [110] and [001] are in the scattering plane. 

The (h/3 k/3 7/2) reflections show identical RXD [2l| 
while the magnitude and sign of the RXD observed at 
the (2/3 -1/3 1/2) reflections depends on I (Fig. |5J). The 
relative magnitude and shape is constant between 10 K 
and 300 K and the diffracted intensity gradually disap- 
pears above the ferroelectric ordering temperature as il- 
lustrated in the inset of Fig. [5] Thomson scattering 
associated with the crystallographic distortion that ac- 
companies the ferroelectric polarization dominates the 
diffracted intensity before and after the edge. The varia- 
tion of XRD with I suggests that the scattering amplitude 
related to the anomalous diffraction of the Fe ions adds 
phase shifted contributions to the scattering amplitude 
from the structural distortion, depending on the I index. 
The structure factor is F = ■ f 3 ■ cxp(iq • Tj) with fj 
the atomic form factor of atom j and r., its position in 
the unit cell. To analyse the energy dependent intensity 
we separate F into an energy independent and an energy 
dependent term, F — F FeLu + F(E)p e , where the first 
term is the Thomson scattering of Fe, Lu and/or O ions 
and the second term is anomalous diffraction due to the 
Fe ions. These terms are written as F Fe Lu Q = A + iB 
and F Fe (E) = ^(f^E) + if'/{E))e lq ^ where A and 
B are the real and imaginary component of the non 
resonant structure factor and f'j(E) and f"(E) are the 
real and imaginary component of the anomalous scatter- 
ing factor of the Fe ions. f'j(E) and f"(E) are related 
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FIG. 2: (a) Intensity of the (1/3 1/3 31/2) reflection as func- 
tion of azimuthal angle 'J', (b) Intensity of the (1/3 1/3 3/2) 
reflection recorded with aa' and an' radiation. The latter 
intensity is multiplied by 10 for clarity. 

through the Kramers-Kronig (KK) relation. 

For the superlattice reflections of this study Fp e {E) 
almost cancels except for the fact that Fe ions in different 
local electronic environments exhibit different energies 
of the Is and the 4p states. This results in different 
transition energies, for example a chemical shift of 4.5 
eV has been reported for f"(E) between Fe 2+ and Fe 3+ 
in water complexes (23j . In addition, ATS arises when 
the Ap states are split. The extended Ap orbitals are 
sensitive to local distortions and orbital order of the 3d 
shell gives rise to splitting of the Ap states, for example 
via the associated Jahn Teller effect. 

Differences in ATS give rise to a modulation in RXD 
intensity as function of azimuthal angle that is related 
to the symmetry of the distortion. Combined with the 
Thomson scattering, whose diffracted intensity is inde- 
pendent of azimuthal angle, this results in a change of 
the relative amplitude. Fig[2^i shows this is not ob- 
served in LuFe204 in contrast to orbitally ordered man- 
ganites [ll|. Furthermore, a significant contribution in 
diffraction with rotated polarization is expected for ATS, 
whereas polarization analysis shows that the diffracted 
intensity with rotated polarization, an' , is weak and ac- 
counted for by the unrotated aa' contribution within the 
polarization resolution of the analyzer (Sec Fig. Wp)- The 
(1/3 1/3 1/2) reflections of Fig. [5] have arbitrary angles 
4>i with the [001] direction (<f> 3 = 73° and 4> 31 = 18°) so 
that ATS is not canceled to zero by symmetry if there 
is orbital order. Moreover, polarization analysis of vari- 
ous other (1/3 1/3 1/2) reflections and several azimuthal 
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angles recorded for 1=29 yielded the same result 

Both aspects signal the absence of local asymmetric 
distortions associated with orbital order and point to a 
rather symmetrical expansion or contraction of the trig- 
onal bipyramids in the ab plane, changing the electronic 
density at the resonant ion while the symmetry of the 
scattering factor is unaltered. 

The energy dependencies of the (ft/3 k/3 1/2) refl.ee- 
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FIG. 3: (a) f"(E) and /2(E) deduced from the refinement of 
RXD (Fig. [5} via the method described in the text and (b) its 
average compared with FDMNES calculation in R3m (dashed 
curve) and (f{'(E) + f%(E))/2 from fi(E) (open circles). 

tions are examined with two different charge states of 
the Fe ions, labeled 1 and 2, and the structure factor at 
(ft/3 k/3 1/2) equals 



Fi{E) 



A' l + iB[ 



+f 1 {E)+if 1 \E)-f 2 {E)-if 2 \E). (1) 



Where A\ and B[ are real constants of arbitrary mag- 
nitude. In particular B[ is non zero due to the charge 
order and breaking of inversion symmetry. The sum 
if" + /2O i s obtained from fi(E) using the optical the- 
orem. In this work we aim to deduce fi(E) and /2(E) 
from the RXD spectra without any assumptions on the 
local distortion. To test the robustness of the specific 
energy dependencies, the series of (2/3 -1/3 1/2) reflec- 
tions are refined with different methods. Besides fully 
independent fi(E) and /2(E), we added the constraint 
f 2 (E - |A) = fi(E + |A), where A equals the chem- 
ical shift. Moreover, we have calculated f"(E) with 
FDMNES using the muffin tin approximation [2l], [24] 



in R3m and used f 2 (E - ±A) = f x {E + ±A) to fit the 
XRD data. The chemical shift between f['(E) and f 2 '(E) 
is similar in the three refinements and equals 4.0(1) eV. 
The first refinement resulted in distinct energy dependen- 
cies for fi(E) and /2(E). This is understood as a result 
of the experimental uncertainty of the XRD spectra and 
minimizing fc{E — ^A) — fi(E + ^A) was added in a 
further refinement to promote similarity between fi(E) 
and /2(E). Figj5] illustrates the resulting fits and f"(E) 
and fi/(E) are presented in Fig. H Both f['(E) and 
/2(E) exhibit characteristic features at similar energies 
above the edge however their magnitudes are different. 
FigH] compares the three different models for the (2/3 
-1/3 7/2) reflection. 

The chemical shift between f"(E) and f%(E) of 4.0(1) 
eV corresponds to the chemical shift between Fe 2+ and 
Fe 3+ in FeO and Fe203 and confirms Fc 2+ /Fc 3+ charge 
order. Besides the invariable chemical shift there are dis- 
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FIG. 4: RXD compared to three models of charge order as 
described in the text. The solid curve corresponds to f" >2 (E) 
as presented in Fig. [3] The dashed and dotted curves cor- 
respond to identical fi(E) and /2(E) except for a chemical 
shift. The dashed curve is fitted to the RXD data whereas 
the dotted corresponds to fi$(E) calculated with FDMNES. 

similarities between the models. The double feature in- 
dicated with B,C in Fig. [5] is not reproduced by the fit 
based on the FDMNES calculations. In addition, the 
broad feature, labeled D, cannot be accounted for with 
f 2 (E - § A) = fi(E + ±A). Yet, features B,C and D 
are either in or out of phase with the Thomson scatter- 
ing (see Fig. [5j. This signifies that the Fe-0 bonds are 
distinct at each Fe site. 

Long range of the Fe 2+ orbitals leads to ATS that is 
distinct from that of the spherical Fe 3+ ion. Such dif- 
ference in ATS between the Fe 2+ and Fe 3+ ions is ab- 
sent and evidence of random orientations of the Fe 2+ 
orbitals. Our findings are consistent with the unconven- 
tional orbital state calculated by Nagano with frustrated 
orbital orientations and large degeneracy in the ground 
state fli ]. 

The orbital state critically influences the superex- 
change interaction and is essential to understand the 
magnetoclcctric properties of LuFe204. In contrast to 
the manganitcs. where the orbital liquid state is discussed 
to be associated with the ferromagnetic and metallic 
state [25j], LuFe2C>4 shows orbital disorder in the ferro- 
electric state. As such it is more appropriate to classify 
LuFe204 as an orbital glass, adding further frustration 
to the already frustrated magnetic interactions of Ising 
spins on a triangular lattice. 

Fe304 and LUFC2O4 both exhibit ferrimagnetic or- 
der combined with magnetoelectric effects [26| but their 
charge and orbital states are evidently distinct. Fe3C>4 
shows merely fractional charge order (l3j and recently 
order of the t 2g orbitals was observed below the Verwey 
transition j2?| . While the superexchange is prevailing in 
Fe304, coulomb interactions dominate in LuFe2 04, re- 
sulting in charge order and crystallographic distortions 
that do not accommodate orbital alignment. The re- 
duced amplitude of /" of Fe 2+ just above the edge (see 
Figl3K) is possibly indicative of this frustration and sub- 
ject of further investigation. 

Finally it is noted that cooling from 380 K in an electric 
field of ~1 MV/m (EFC) did not result in a significant 



change of the (/i/3 k/3 1/2) superlatticc reflections which 
contradicts the suggestion that ferroelectric charge order 
is the ground state only after EFC [H]. 

In conclusion, our RXD data show an almost complete 
Fe 2+ /Fe 3+ charge order as origin of the superlattice in 
LuFe2C>4. The Fe scattering factors are isotropic and 
in agreement with frustrated and random orientations of 
the Fe 2+ orbitals, forming an orbital glass state. 
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Supplement 

In this supplement we present additional details of our 
study that might be of interest for the specialist reader. 
In the first section we present details of the FDMNES cal- 
culations, in the second section details of resonant x-ray 
diffraction (RXD) recorded at other superlattice reflec- 
tions and in the third section details of the experimental 
geometry. 

1. FDMNES 

The FDMNES calculation presented in Fig 3. a of the 
paper was performed using the muffin tin approximation 
[l| with space group Rim (no. 166), crystal parameters 
a = 3.4406, b = 3.4406, c = 25.280, a = 90°, (3 = 90°, 7 
= 120° and atom positions of Lu at (0,0,0), Fe at (0, 0, 
0.21518) and O at (0, 0, 0.1291) and (0, 0, 0.2926). The 
result converged at a cluster radius of 6A. 



1 

- LuFe 2 4 
- 


1 1 1 

(1/3 1/3 7/2) 

« (2/3 -1/3 7/2) 

| - (1/3-2/3 7/2) 
\ . (-1/3-1/3 7/2) " 
1 f; - (-1/3 2/3 7/2) 
J H (-2/3 1/3 7/2) 


\ A 1 

- vyi 

u 

V 

1 


i i 1 



7100 7120 7140 7160 7180 
Energy (eV) 

FIG. 5: Normalized integrated resonant diffraction intensity 
of the charge order (ft/3 fc/3 7/2) reflections recorded with 2D 
detector at T=15 K, corrected for background, polarization 
and absorption. 



2. IN PLANE SYMMETRY 

The symmetry of the ferroelectric structure of LuFe2C>4 
is not unambiguously known at this stage and there are 
possibly several tensors of x-ray susceptibility (ATS) al- 
lowed. These terms are significant in systems where 
charge order is connected with orbital order, such as man- 
ganites and magnetite, resulting in strong azimuthal an- 
gle dependences as well as observable contributions in the 
rotated channels [1, 01]. In contrast, charge order without 
orbital order as observed in for example NdNi03,results 
in RXD that is independent of azimuthal angle [4|. The 
(1/3 1/3 31/2) and (1/3 1/3 3/2) reflections have angles 
(j) = 18° and 73° with the [001] direction so that ATS 
is not canceled to zero by symmetry if there is orbital 
order. Moreover we recorded azimuthal angles of 0°, 30° 
and 60° for the (1/3 1/3 29/2) reflection (</> = 19°) and 
polarization dependencies for various (1/3 1/3 1/2) reflec- 
tions with Z=23-37 (0 = 23°-15°) (at beamlinc ID20 of 
the ESRF) which yielded the same result. 

Figure \5\ illustrates resonant X-ray diffraction (RXD) 
for (h/3 k/3 7/2) reflections (0 = 54°) recorded at the MS 
beamline of the SLS. The RXD of these six reflections is 
identical, consistent with eq. (1) in the paper. 

3. SURFACE DIFFRACTOMETER 

RXD as illustrated in Fig. 1 of the paper and Fig. [5] 
of this supplement were recorded at the surface diffrac- 
tometer of the X04SA beamline Q using the Pilatus 2D 
detector Q. This station is a large 2+3-circle surface 
diffractomctcr with two circles for the sample and three 
for the detector, plus a hexapod for precise alignment 



of the sample surface relative to the diffractometer axes. 
The diffractometer was operated with a vertical sample 
surface orientation. 

Horizontal polarization of the incoming radiation was 
used which corresponds to n incident polarization when 
the scattering plane is horizontal. When the scattering 
plane is not horizontal the incident beam contains both 
a and n polarized components. The polarization of the 
incident radiation is most conveniently characterized by 
the Poincare vector Pf which equals +1(-1) for a (it) 
incident radiation 0. Table U gives P^ for the super- 
lattice reflections presented in this study. Table Q] also 
summarizes the angles that define the orientation of the 
sample (a, oj v ) and detector (7, 5) with respect to the 



TABLE I: Orientation of the sample (a, uj v ), detector (7, 
S) and the polarization of the incident radiation P^ for se- 
lected superlattice reflections recorded at the surface diffrac- 
tion endstation of beamline X04SA with incident x-ray energy 
7112.4 eV. 
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Pc 


(1/3, 1/3, 7/2) 


10.00 


66.29 


14.09 


19.73 


0.37 


(1/3, -2/3, 7/2) 


10.00 


-48.32 


12.97 


20.31 


0.46 


(-1/3, -1/3, 7/2) 


10.00 


-106.19 


14.39 


19.29 


0.33 


(-1/3, 2/3, 7/2) 


10.00 


127.84 


15.70 


18.44 


0.21 


(-2/3, -1/3, 7/2) 


10.00 


-169.65 


15.48 


18.27 


0.21 


(2/3, -1/3, 7/2) 


10.00 


10.06 


12.69 


20.49 


0.49 


(2/3, -1/3, 11/2) 


10.00 


11.39 


21.03 


20.99 


0.07 


(2/3, -1/3, 13/2) 


10.00 


14.56 


25.27 


21.05 


-0.10 


(2/3, -1/3, 15/2) 


10.00 


17.15 


29.75 


21.43 


-0.23 


(2/3, -1/3, 17/2) 


10.00 


20.69 


34.61 


20.83 


-0.38 


(2/3, -1/3, 21/2) 


10.00 


30.45 


43.68 


20.01 


-0.56 
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incoming radiation. The sample is mounted at the centre 
of rotation and a and 7 are the angles in the horizontal 
plane between the incoming beam and the sample sur- 
face and detector position respectively. & v is the sample 
rotation about an axis normal to the sample surface and 
(5 is the angle between the horizontal plane and the de- 
tector position respectively. No polarization analysis of 
the diffracted beam was performed and the pixel detector 
records both rotated and unrotated polarizations. 
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